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Abstract 

Liquid crystals have attracted enormous interest because of the variety of their phases and 
richness of their application. The interplay of general physical symmetries and specific molecular 
features generates a myriad of different phenomena. A surprising behavior of liquid crystals is the 
reentrancy of phases as temperature, pressure, or concentration are varied. Here, we review the 
main experimental facts and the different theoretical scenarios that have guided the understanding 
of bulk reentrant nematics. Recently, some computer simulations of a system confined to nanoscopic 
scales have found new dynamical features of the reentrant nematic phase. We discuss this prediction 
in relation with the available experimental evidence on reentrant nematics and with the dynamics 
of liquids in strongly confined environments. 
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I. INTRODUCTION 



In 1888 the Austrian chemist and botanist Friedrich Reinitzer observed that the organic 
compound cholesteryl benzoate exhibits two melting points. He found that it first melts 
at 145.5 °C to a cloudy liquid, and then at 178.5 °C it becomes suddenly c .ear. Later, the 
pioneering studies of Lehmann [1], Schenk j2|, Vorlander jsl and Friedel 4| recognized a 
new state of matter intermediate between liquid and solid. Liquid crystals (LC's) were 
discovered. We now know a multitude of different liquid crystalline phases, roughly divided 
into two families, (i) Thermotropic LC's undergo phase transitions as the temperature 
(T) is changed, (ii) Lyotropic LC's exhibit different phases as the solvent concentration is 
changed. The most commonly encountered phases are nematic (N), where the system has 
a preferential orientation given by the nematic director n, but the centers of mass have no 
long-range correlation (see the bottom right panel of Figil]); smectic (S), where molecules 
organize in layers (see the bottom left panel of Fig{T]); columnar, where two-dimensional 
order is established (see Fig|2]); cholesteric, where the average orientation spirals with a fixed 
pitch as different planes are traversed, and blue phases, which consist of three-dimensional 
arrangements of cholesteric tubes. 

Many different molecules form liquid crystalline phases, and their modern importance is 
difficult to exaggerate. Liquid crystal displays for electronic devices are ubiquitous. Because 
of the anisotropy of LC molecules and of their interactions, LC's can easily polarize light 
passing through them. Highly versatile lasers can be made with LC systems, promising 
numerous applications jsj. Liquid crystals are extremely important also for biology. Some 
examples (far from being comprehensive) are the following: The lipid bilayer of a cell mem- 
brane is a stable LC phase on which depends the very existence of cells; strands of DNA 



are also known to form LC's 



,|7j; spider's silk has exceptional materials properties due to 
the parfcular LC organisation of its proteins H; finally. LC's are ushernrg a technological 

revolution with multiple applications in biology and medicine as biosensors [9j. Many ex- 
tensive reviews of the state of the art of LC crystals exist, see for example [lo| for modern 
applications of LC in self-assembling supramolecular structures. 

Among the interesting features of LC's that make them so versatile there is a surprising 
sequence of phase transitions of rod-like molecules that still poses challenges to both experi- 
mental and theoretical studies. The typical sequence of LC phases encountered starts at high 
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T with an isotropic (I) phase, where the system is characterized by complete translational 
and orientational invariance. As T decreases a preferential molecular orientation emerges - 
the nematic director - breaking the rotational symmetry. At even lower T, S layers appear 
indicating that now translational invariance is broken, too. 

Surprisingly, in 1975 Cladis discovered [llj] that when two compounds, whose molecules 
have two benzene rings and a strongly polar cyano group, were mixed together, a N phase 
formed at T above the S phase (which is usual), but also at T below the S phase (which is 
surprising). The N phase at lower T was termed reentrant nematic (RN). The presence of a 
RN phase is surprising because in the vast majority of cases as T is lowered the symmetry 
of the stable thermodynamic phase decreases, i.e. the state of the system is invariant under 
a decreasing number of symmetry operations, eventually becoming a crystalline solid. 

After the seminal work of Cladis ll|], other systems in a variety of different conditions 



were found to have a RN phase. Today, RN behavior is a well-established phenomenon. 
Here, we focus on the past and present understanding of RN, with special attention to the 
dynamics of RN phases. 

Naturally, nearly 40 years after Cladis' discovery of the RN phase one may wonder: Is 
there still something to be understood about these phases? It is the primary purpose of 
this article to demonstrate that it is indeed so. For exmaple, recent numerical investigations 
on the dynamics of LC's confined in nanoscopic pores suggest some interesting dynamics 
features. In particular, it was found 12|, |l3| that diffusivity in the direction of the molecular 
long axis is greatly enhanced. This appears to be an analogous situation to the levitation ef- 



'ect in zeolites and carbon nanotubes 



14 



15| . There is some favorable experimental evidence 
16| supporting the theoretical predictions of enhanced parallel diffusivity 12, ist. However, 



the question is not yet settled. Our critical review of the available experimental work on RN 
dynamics reveals that numerous controversies are not yet resolved. A main problem identi- 
fied here is that reentrancy is not caused by a single microscopic mechanism, but rather by 
many. Their identification is still an open question. Consequently, the investigations of the 
dynamical behavior of RN's are often divergent. 

The remainder of this Review is organized as follows: in Sec. |TT] we review the seminal 
experiments that discovered RN phases in both polar and nonpolar substances. In Sec. IIIII 
we discuss the main experimental investigations on the nature and properties of RN's, while 
in Sec. IIVI we review the different theories and scenarios put forward to rationalize RN 
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FIG. 1: Two examples of mesophases formed by rod-like molecules. Smectic A (bottom left) and 
nematic (bottom right). 

existence and features. In Sec. |V]we review the more recent computational studies of RN's, 
and, finally, Sec. |VT] lists our conclusions. 



II. HISTORY 



The first experimental detection of RN in a LC is reported by Cladis There, an 

amount of hexyloxybenzilidene amino benzonitrile (HBAB) was mixed with cyanobenzili- 
dene butylaniline (CBOOA). Individually, the former has a N phase for 35 °C < T < 101 °C, 
and the latter has a Sa phase for 44 °C < T < 83 °C and a N phase for 83 °C < T < 108 °C. 
When mixed together, as T is lowered a sequence of I, N, Sa, and RN is found for a range of 
molar fractions c; the curve bounding the Sa phase in the T — c plane can be approximately 
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FIG. 2: Two examples of mesophases formed by disk-like molecules. Columnar (bottom left) and 
nematic discotic (bottom right). 



fitted with a parabola 

c = co + /3(TNs-To)^ (1) 

where cq ^ 0.09, /3 ^ -5.4 x 10^ °C~^ Tq ^ 61 °C, and Tns is the N-Sa transition T 
measurements of the bend elastic constant ll| constitutes the first examination of possible 



The 



structural differences between the N and RN phases. In this early experiment it was shown 
that the elastic constant has a very similar T-behavior as the Sa phase is approached from 
above or below pLL]. Cladis emphasized that CBOOA has an incommensurate Sa phase, 
that is, the spacing between S layers is not an integral multiple of the molecular length. A 
possible way to understand how an incommensurate spacing can arise is to imagine a S layer 
composed of dimers. Because of the strong polar character of cyano compounds molecules 
naturally tend to form dimers to minimize the interaction energy; these dimers in turn 
organize themselves in S layers, or, in other words, each S layer is in reality a double layer of 
molecules. This initial observation is consistent with other experiments [l7[ and has proven 
itself fruitful in the understanding of the structure of RN phases. 

Following the assumption of the double-layer nature of the Sa phase a series of pure 
compounds and some mixtures were studied [l8|. Upon increasing pressure, RN phases were 
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found in the supercooled region of polar substances such as 4-cyano-4'-octyloxy biphenyl 
(COOB) or CBOOA. What these substances have in common is the amphiphilic nature 
of their molecules. Unlike most LC's the polar part (the aromatic rings) is not in the 
middle of the chains, but rather at one end. The polar cyano group has a dipole moment 
of 4.5 D. The nonpolar part (aliphatic chain) forms the opposite end of the amphiphile. 
This structure will naturally favor a certain degree of dimerization, since the polar groups 
will preferentially interact with each other through long-range forces, and the nonpolar tails 
through short-range forces. Cladis et al. [l8| proposed that it is the short-range interaction 
between nonpolar tails what stabilizes the S layers. This hypothesis is successfully tested by 
measuring the dependence of the maximum stability pressure Pm of the Sa on the number of 



carbon atoms in the nonpolar tails. It is 
methylene groups increases. It is argued 



bund that Pm increases linearly as the number of 
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19| that as pressure increases (or equivalently 



as T decreases) the interaction between the polar groups becomes repulsive, and at the same 
time the nonpolar tails are somewhat compressed, lowering the energy barrier to permeation 
through the S layers. Both effects lead to a destabilization of the S order, and hence to the 
formation of a RN. 

After 1979, the occurrence of RN could not be dismissed as a peculiarity of few molecules, 
since a stable RN at atmospheric pressure was found for molecules exhibiting three benzene 



rings, and, importantly, not in mixtures but in^pure compounds 
reentrant smectic phases were also found 



23 



20|-|22j. Additionally, 



24| . In the following years different groups 



found many different systems exhibiting RN phases (see, e. g., t he work of Sigaud et al. 22 1 



and references therein), and also multiple reentrant phases 25 



, l26| . It has been speculated 



that quadruple reentrance is possible if compounds with four or five benzene rings are 



employed 



271]. Cladis's own account of the discovery and thermodynamic understanding of 



RN phases provides many more details 



28|. 



Until 1983 it was believed that only molecules with a strong polar group showed a reen- 
trant behavior. However, Diele et al. showed the existence of RN phases in terminal- 



nonpolar substances 29|, |30|. Importantly, their X-ray measurements showed that the thick- 



ness of the S layers is equal to the average molecular length. Hence, dimerization cannot take 
place in these nonpolar compounds and the occurrence of a RN phase must have a different 



physical origin 



29, 



30| (see Sec. |IV] for the theories applicable to nonpolar systems) 



A question naturally arising is whether RN and N phases are different. One can consider 
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two aspects: structure and dynamics. In the following we explore the experimental evidence 
addressing this question. We can, however, anticipate that, since pure compounds and 
mixtures, polar and nonpolar molecules all exhibit RN phases, the picture emerging will be 
rather complex and diverse. 



III. EXPERIMENTAL INVESTIGATIONS 



s 

ing 



;ructural differences of the N and RN phases were investigated with X-ray scatter- 



31| for mixtures of N-(]9-hexyloxybenzylidene)-^)-aminobenzonitrile (HBAB) and N-p- 



cyanobenzylidene-j5-n-octyloxyaniline. No qualitative difference was found between the N 
and RN phases, i.e. the Bragg peak indicating the layer thickness appears always in the 
range 33 — 35A; however, for some mixtures the RN was found to coexist with microcrystal- 
lites. Guillon et al. concluded that these crystalline fluctuations share the same structural 
properties and are energetically similar to the RN phase j3l| . 

Nuclear magnetic resonance (NMR) offers a number of techniques to investigate both 
structure and dynamics. The first NMR study of a pure LC system with a RN phase 



was carried out by Miyajima et al. 



32| . For pure 4-(4"-octyloxybenzoloxy)-benzylidene-4'- 



cyanoaniline (OBBC) they found a smooth continuous increase in the orientational order 
parameter across the N, Sa, and RN phases. However, when studying a binary mixture with 
a spin-probe impurity probe they found a steep increase of the nematic order parameter in 
the RN phase. While proton NMR is sensitive to the aromatic core, the impurity probe 
method is not. Therefore, the conflicting results in the orientational order parameter can 
be ascribed to the different techniques. From the T behavior of the spin-lattice relaxation 
time Ti measured at 29.8 MHz for OBBC Miyajima et al. conclude that translational 
self- diffusion plays a dominant role in Ti relaxation in the RN phase (with a characteristic 
Arrhenius T- dependence), while the N phase Ti behavior is dominated by nematic director 
fluctuations. 

However, for a LC mixture Dong reached the opposite conclusion that self- diffusion ap- 
pears to be the relevant relaxation mechanism in the N phase, while nematic director fluc- 
tuations dominate the RN phase^ 33 -|35[| . Using proton NMR he studied the T dependence 



of Ti for a LC binary mixture 



34| and found an Arrhenius behavior in the N phase. 



and no T-dependence in the RN phase within experimental error; that director fluctuations 



7 



dominate the relaxation of the RN phase is further corroborated by the hnear relationship 
between T^^ and where to is the angular frequency, as predicted by the theory of spin 

relaxation through a director-fluctuation mechanism 35 1. 

We cannot avoid stressing at this point something that has not been clearly stated in 
previous works. Although reentrancy is found in many different substances, their dynamics 
need not be the same. Specifically, pure compounds and mixtures do not necessarily have 
to exhibit the same dynamical properties. Even though there are discrepancies in the ex- 
perimental results, we can safely conclude that the dynamics of N and RN are markedly 
different. 



A somewhat intermediate situation is described by Sebastiao et al. [36[ who performed 
proton NMR experiments on a polar LC compound which exhibits N, partial bilayer S, RN, 
and reentrant Sa phases. Self-diffusion is shown to contribute little to Ti in the N and 
RN phase. Instead, fluctuations of the director and reorientations are the main dynamical 
processes contributing to Ti. Interestingly, however, the ratio between the translational 
self- diffusion parallel and perpendicular to the director is slightly higher in the RN phase 



than in the N phase [36[. Those results [36| are consistent with a dynamical process of 
association and dissociation of molecular groups, which are important in the double-layer S 
and N phases. 

By using ^^^Xe NMR Bharatan and Bowers [stI] studied two LC mixtures and found that 
Ti (and also the spin-spin relaxation time T2) shows Arrhenius behavior in both N and RN 
phases. The activation energies extracted from the data do not behave in a similar fashion 
for the two systems studied. For a binary mixture, the activation energy in the RN is more 
than two times the value in the N phase, whereas for a ternary mixture, there is almost 
no difference in the activation energies. The increase in activation energy for the first LC 
mixture is ascribed to changes in the molecular packing in the RN phase 37l |. 

Although a number of other studies were carried out on mixtures and pure compounds 
with a RN phase 38-4^ using deuteron NMR little can be said about translational self- 
diffusion, because this technique is sensitive to other types of motion, e.g. intramolecular 
rotations, reorient at ional dynamics, or supramolecular motion, such as order director fluc- 
tuations. The picture arising from the work of Dong and collaborators is one of very subtle 
changes between the phases, and especially between N and RN phases. 

The structure of RN's also poses some unanswered questions. For example, is there a 
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change in nematic order parameter Q at the transition from S to RN? The work of Vaz et 



al. 



411 ] shows that Q increases upon entering the RN. They measured the T-dependence of 



Q in a LC binary mixture with deuterium NMR and found that at the Sa-RN transition 



Q displays a sharp increase in slope 41j. To rationalize this behavior they considered an 
expansion of the free energy difference AF = F/v + Fs — Fi between the N, S and I phase 
in terms of Q and a smectic order parameter \E', that to leading order is 



(2) 



where a = ao(T — Tjn), a = ao(T — Tna), Tjn and T^a are the I-N and the N to Sa 
transition T, respectively, and all the constants are positive except /3 < 0. This expansion 
of AF produces a first order I to N phase transition, and a second order N to Sa transition. 
Following Cladis 42|, they assumed a quadratic coupling of the order parameters propor- 
tional to Q^'^^ to describe the RN phase and derived that upon entering the RN from the 
Sa phase the nematic order parameter should be positively perturbed ( 6Q > 0) such that 
the phenomenological description is in agreement with the experiment s [41 1. 

However, using deuterium NMR on a binary mixture, Emsley et al. |43| found only subtle 
changes in the order parameters at the Sa-RN phase transition. Contrasting results are also 
presented by Dong et al. 4^] , where deuterium NMR experiments showed no enhanced ori- 
entational order in the RN phase, in contradiction with both a Landau free energy expansion 
theory and a McMillan theory of RN's. 

Electron spin resonance experiments on the LC mixture of hexyloxy biphenyl (60CB) 
and octyloxy biphenyl (80CB) found no dramatic structural changes at the Sa-RN transi- 
tion using a number of different spin probe molecules and similar dynamics of the probes in 



reentrant and nonreentrant LC 
et al. 



45| . Nayeem and Freed also concluded, contrary to Dong 



that their results are "not consistent with saturation in pair formation being a 
necessary precursor to reentrance" j45|. The differences in these experimental results can 
be ascribed to different experimental techniques sensitive to a different extent to the orien- 
tational degrees of freedom of the aromatic cores or of the molecular chains. Nonetheless, 
the question still lingers: What is the behavior of the nematic order across the S-RN phase 
transition? 

A very different experimental approach was adopted by Ratna et al. [l6^ who measured 
the electrical conductivity in the direction parallel (ctn) and perpendicular {a±) to the optic 



axis, which in the orientationally ordered phases coincides with the nematic director. They 
studied the pure compound 4-cyanophenyl-3'-methyl-4(4'-n-dodecylbenzoyloxy)benzoate (12 
CPMBB) and the mixture 60CB/80CB. For the single component system they found a 
remarkable increase of the ratio u\\/iy_i_ as the system approaches the RN-Sa transition tem- 
perature TjiN-A- At the lowest T studied for the pure compound (1.2 °C below Tji^^a) o'w 
is approximately equal to 16cr_L. For the mixture 60CB/80CB, they found the same quali- 
tative increase of the ratio cr\\/a±, but at the lowest T studied in the RN phase the value of 
this ratio was only 1.8 [16]. Although the results for the mixture would require additional 
analysis, what was found for the pure compound is strong evidence supporting a scenario of 
increased mass transport in the RN phase. 

Another insight into the differences between N and RN phases comes from the work of 



Nozaki et al. 



46| who found that for the mixture 60CB/80CB the rotations around the 



short molecular axis are greatly hindered in the RN phase with respect to the N phase. 
Specifically, by measuring the T-dependence of the complex permittivity with time-domain 
reflectometry, they found a higher value of the activation energy in the RN compared with 
the N phase. 

There are limited modern investigations of the RN state. Yethiraj et al. measured proton 
NMR spectra of solutes dispersed in a mixture of 60CB/80CB to determine the smectic 
and nematic order parameters in the Kobayashi-McMillan theory [47 1. Das and Prasad 48] 
measured the rotational viscosity of a LC binary mixture and found similar values in the N 
and RN phases. However, there are strongly contrasting results 0, Q suggesting that the 
rotational viscosity of RN may be two or three orders of magnitude higher than in the N 
phase. 

Finally, hydrodynamic flow in RN's has received some attention. Bhattacharya and 
Letcher measured capillary shear flow for a ternary mixture with a RN phase. They found 
that the N and RN phases have identical flow properties, and the hydrodynamic theory of 
N phases can be applied to RN phases [49|. Also, Jadzyn and Czechowski [50| measured 
the viscosity for the well studied mixture 60CB and 80CB, and found a distinct decrease 
of viscosity upon entering the RN phase. 
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IV. THEORETICAL SCENARIOS FOR REENTRANCY 



The experimental findings discussed in Sec. UTIlwell illustrate the complexity of RN phases. 
It is therefore not surprising that a similar rich picture also emerges from theoretical stud- 
ies. Prost and Barois derived a Landau theory of LC exhibiting reentrancy from general 



considerations of symmetry and on the basis of the relevant length scales involved |5ll . |52 |. 
In a S phase a density modulation appears. On average the molecular centers lie on parallel 
planes. Thus, it is natural to consider the expansion of the density in a Fourier series 

p = Po + ^(#^''-^) + ... (3) 

where the wave- vector q = (27r//)n, / is the layer thickness, and n is the director. However, 
polar substances and mixtures exhibiting RN phases show two incommensurate periodicities 
in X-ray scattering experiments, i.e. one periodicity equal to the molecular length /, and 
a second one /' that is not an integer multiple of the first, typically /' ~ 1.2 - 1.5/. It is 
necessary to introduce the total dipole moment P(r) to describe a system of polar molecules 

P(r)^l5^M(r-r,), (4) 

i 

where pi is the dipole moment of molecule i. A potential can be derived from P(r) that 
may be used as a second order parameter 

P(r) = i-V<^. (5) 

With these two order parameters, p and </>, a Landau- Ginzburg free energy functional can 
be written as 

F = 1 1 dh{A,p' + A^<i? - Dp<p' + li?,p4 + + Cp^^} . (6) 

From the competition between the different length scales contained in Eq. ([6]), Prost and 



collaborators 



51 



52| obtained a phase diagram with different S and also RN phases. 
The first microscopic theory of the physical origin of RN phases was elaborated by Berker 
and Walker in 1981 [53]. They started from the observation that most reentrant nematogens 
are composed of molecules with a strong polar group at one end and an attached aliphatic 
chain forming the other end of the molecule. The dipolar interaction is assumed to be 
antiferroelectric. Considering a plane normal to the nematic director, antiferroelectric long- 
range order cannot be established in a periodic regular arrangement of molecules. However, 
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Berker and Walker argued that in a liquid molecules can relax their positional degrees 
of freedom and configurations of triplets with two short antiferroelectric bonds and one 
unfavorable ferroelectric, but longer bond will be possible. In this way frustration is avoided 
and the system can condense to a liquid of dimers forming a two-dimensional network. Hence, 
a bilayer S phase forms. As T is lowered positional order arises that leads to increasing 
frustration and to the breaking of the dimers. As a consequence a nematic phase appears 
again. The interaction of two dipoles Si and S2 separated by a distance |ri2| is described by 
the potential 



ri2 



where |sj| = 1. For B < A the interaction is antiferroelectric, whereas for B = A it is 
purely dipolar, and for B > A the interaction becomes ferroelectric. By integrating out 
the positional degrees of freedom and distinguishing different degrees of interaction, Berker 
and Walker mapped the initial system onto an Ising model with annealed bond disorder. 
The phase diagram obtained |53| shows a RN phase for a range of model parameters and 



also a doubly reentrant phase sequence. Berker et al. 



27 



54j extended their theory to 



describe multiple reentrant behavior, considering partial bilayer (SAd) and monolayer (Sai) 
smectic phases, i.e. N-Sah-N-Sah^N-Sai , which was found experimentally in at least one 
pure substance and its mixture [55| and Sc phases [56|. The scenario of frustrated dipolar 



interaction 



27 



54| (commonly referred to as "spin-gas theory") has played a major role in 



understanding different instances of reentrancy in LC 



28|. An implication of the scenario 



described above is that the N-Sa transition is different from the Sa-RN transition. This is 
at odds with some experiments that show very subtle or no structural difference between 



the nematic phases for LC mixtures [45|, |57| and a pure compound 



A different approach to a microscopic theory for reentrancy in nematogens was proposed 



by Longa and de Jeu 



59|: Pairs of molecules with antiparallel dipoles associate to form 



dimers in dynamical equilibrium with monomers. As T is decreased the fraction of dimers 
in the LC system increases and dipolar dispersive forces between the dimers lead to the 
formation of a Sa phase. It is important to note that monomers under these conditions 
help stabilizing the Sa phase because monomers can fill the voids left by the bulkier dimers. 
However, as the dimer fraction increases it becomes increasingly difficult to pack them with- 
out perturbing the S layers. It is instead, as Longa and de Jeu argue jsol, entropically more 
favorable to disrupt the S layering by forming a RN phase. By extending McMillan's theory 
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of the N-Sa transition [60], Longa and de Jeu used a mean-field approach to solve for the 



one-partic 
rators 27, 



e distribution functions of monomers and dimers. Although Berker and coUabo- 



53 



54 



56| too considered dipolar interactions to be fundamental, they viewed the 



S phase as a percolating network of antiparallel pairs where frustration is avoided only be- 
cause of liquid positional disorder. In our view, contrary to the conclusions of Figueirinhas et 
al. 171], it is difficult to exclude either theory in favor of the other, because the fundamental 



experimental predictions of the spin-gas theory, and of the McMillan-like theory are quite 
similar: They both predict a roughly T-independent value for the smectic layer thickness 
and antiparallel associated pairs. 

It is worth mentioning other approaches to explain reentrancy in LC. For example. Luck- 



hurst and Timini 



6l| found a N — )■ Sa — )■ RN phase sequence upon decreasing T. They only 



assumed a linear T dependence of the parameter a that in the original McMillan theory [60 1 
is empirically related to the length of the molecular alkyl chain. In practice, as a decreases 
with decreasing^ T, the particles described by the theory turn from dimers to monomers. Fer- 
rarini et al. [62] showed that reentrant phases can also be recovered by considering chemical 
reactions of isomerization and dimer formation under reasonable conditions. 

All the models and theories discussed above refer to polar molecules. It it natural that 
most theoretical work has focused on this class of LC's, because experimentally reentrancy is 
predominantly found in LC molecules with strong dipolar groups (typically cyano). However, 
RN phases have also been observed in nonpolar substances 29|, l30| . Dowell |63| devised a 
lattice theory for a single component system that focuses on molecular chain flexibility. In 
Dowell's model the S phase is formed because of segregated packing of cores and tails. As 
T is lowered, however, the molecular chains become increasingly more stiff, destabilizing 
the S layers. Eventually, it becomes entropically more favorable to disrupt S layering such 
that a RN phase is formed. X-ray [57] and ESR 45| experiments support Dowell's scenario 
of reentrancy in systems that do not show signs of dimerization 29|, |30[ • Bose et al. 
extended McMillan's theory to include chain flexibility and also obtained a RN for a single 
component nonpolar system. 
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V. COMPUTER SIMULATIONS 



In modem times vast possibilities of theoretical investigation have been ushered by the 
growth of computational power and the sophistication of simulation algorithms. More and 
more complex fluids and microscopic structures can now be simulated 65| . Computer simu- 
lations were first applied to the study of RN's by Netz and Berker 56(|. They solved numer- 
ically the lattice spin-gas model introduced in Sec. llVl with Monte Carlo (MC) calculations. 
The first off-lattice MC simulations were performed only in 2005 by de Miguel and Martin 



del Ri'o 



66| . They considered a single-site molecular model for mesogens, where molecules 



are represented by parallel hard-core ellipsoids with a spherically symmetric square well 
potential (around each molecular long axis). In that study a RN phase and two tricritical 
points are found, one located on the N-Sa transition line and the second one on the Sa-RN 
transition line. In this simple model, the presence of a RN in a nonpolar LC is rationalized 
in terms of entropy. In the Sa phase the energy is minimized and within the fluid state no 
further decrease is possible. However, the free energy can still be minimized by an increase 
in entropy. Thus, the loss of positional order is entropically driven 66j. 

Though the system of parallel hard-core ellipsoids [66| offers a clear picture of the origin of 
reentrancy in nematics, it has the shortcoming of neglecting rotational dynamics altogether, 
because the molecular long axes point in the same fixed direction. Suppressing orientational 
dynamics leaves open the question of whether the RN in such simple models is robust against 
orientational fluctuations. 



A. Confined Mesogens 



Confined fluids have become an extremely active field of research in recent years 67H69|. 
Novel phenomena take place in liquids confined to narrow spaces, e.g. wetting, layering 
and many others induced by the interaction with the confining surface. The interplay of 
these new phenomena with the physics of bulk fluids gives the possibility of both testing 
our understanding of molecular mechanisms, and devise concrete applications not otherwise 
feasible. In fact, applications of confined fluids range from nanotechnology to biomedical 
devices. It is therefore natural that LC properties and phase transitions come under a 
new scrutiny. Furthermore, nano confinement is a physically sensible (and realistic) way of 
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inducing ordering fields into the fluid system. 

Below, we review computer simulations of a pure LC system confined between walls 
separated by a nanoscopic distance (few tens of molecular diameters) performed by the 



present authors and their colleagues 



12 



13| . These are the first simulations to fully take 



into account a three dimensional fiuid with rotational degrees of freedom. 



1. Model 



Because the focus is on the dynamical behavior of the LC's, molecular dynamics (MD) 
was employed. Specifically, we Q, Q adopted the Gay-Berne-Kihara (GBK) model for 
prolate mesogens 70|. The GBK model potential conveniently takes into account both the 
anisotropy of mesogen interaction and the spherocylindrical molecular shape 7l|] which is 
considered to approximate prolate mesogens better than ellipsoids 70 1. 

In the GBK model the interaction between a pair of spherocylinders i and j depends on 
the molecular orientations represented by the unit vectors ttj and itj, respectively and their 
distance Vij = Vi — Vj. Specifically 

12 



Us = 4:es{rii,Ui,Uj^ 



a 



a 



(8) 



where r = r/r, r = \r\, and the function dfj{rij,Ui,Uj) is the minimum distance between 
that pair of molecules j72|. The orient at ion- dependent interaction strength in Eq. (|8]) is 
given by 

\2 



es[rij, Ui, Uj) 



2 



1 + x'Ui ■ Uj 



yfij ■ Ui 



Vij ■ Uj 



\2n 



X'Ui ■ u 



X 



X^{u, ■ Ujf 



(9) 



where the parameters x x' ^^e given by 



X = 



1 



X' 



/«2 + 1 



(10a) 
(10b) 



1 ■ 

In these last two expressions parameters k = L + 1 (L in units of cr), and k' represents 
the interaction strength for a side-side relative to an end-end configuration of a pair of 
sphero cylinders . 
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We consider a LC system confined along the z direction by two planar, structureless walls. 
We model the fluid-substrate interaction via 



10 / \ 4 

a \ a 



(11) 



where the parameter efs is the strength of fluid-substrate interaction and ps = 2/£^ is the areal 
density of the substrate where i/a = 1/ is the lattice constant of a single layer of atoms 
arranged according to the (100) plane of the face-centered cubic lattice [l2|. The diameter 
a of these substrate atoms is taken to be the same as the diameter of a sphero cylinder of 
the conflned fluid phase. 

As already mentioned, a solid substrate is also a physically realistic way of inducing an 
ordering field in a LC system. A number of experimental techniques were devised to realize 
specific substrate features, and especially anchoring. Some example are rubbing, polish- 



ing, exposing to photolithographic relief, or evaporation of oxide films [73|. In Eq. flTT]) . 
< < 1 is the so-called "anchoring function". It allows one to discriminate energeti- 
cally different orientations of a molecule with respect to the substrate plane. We employed 
degenerate planar anchoring [tJ] 

g{u) = {u ■ + {u . (12) 

where (a = x, y, z) is a unit vector pointing along the a-axis of a space-fixed Cartesian 
coordinate system. Hence, any molecular arrangement parallel with the substrate plane is 
energetically favored, whereas a homeotropic alignment of a molecule (n || e^) receives an 
energy penalty by "switching off" the ffuid-substrate attraction altogether. 

The MD simulations were carried out using a velocity- Verlet algorithm for linear molecules 
After equilibrating the system to the chosen pressure and T, the fluid system is sim- 



ulated in the microcanonical ensemble to obtain information on its dynamics [13(]. We 
simulated systems containing up to = 1500 molecules. Quantities of interest will be ex- 
pressed in the customary dimensionless (i.e., "reduced") units. For example, length will be 
expressed in units of a, energy in units of eg, temperature in units of eff/fce, time in units 
of (cr^m/eff)"'^/^ using m = 1, and pressure Py in units of k^T /a^ where P\\ = |(Pxx + Pyy) is 
related to diagonal components of the pressure tensor P acting in the x-y plane. 
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FIG. 3: Color online) Plots of nematic S* (•) and smectic order parameter A (■) as functions of 
pressure Py at (a) T = 4.0, and (b) at T = 6.0. 



2. Results 



To characterize the structure of N, Sa, and RN phases a suitable order parameter is the 



so-called alignment tensor defined as 



76 



1 ^ 

Q = (13) 

i=l 

where "®" denotes the direct (i.e., dyadic) product and 1 is the unit tensor. Hence, Q is 
a real, symmetric, and traceless second-rank tensor which can be represented by a 3 x 3 
matrix. Its largest eigenvalue is commonly used as the N order parameter S. 
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To characterize the layering of the fluid in the direction of n characteristic of S phases a 
suitable order parameter is given by the leading coefficient of the Fourier series of the local 
density p{r) 



A 



d 



(14) 



X^exp 

1=1 

where d is the spacing between adjacent Sa layers. From Eq. ( [T^ it is also apparent that 
A G [0, 1], where A 0, and A ^ 1 in ideal N and Sa phases, respectively. 

Plots of the dependence of S and A on P\\ are presented in Fig. [3] for T = 4.0 and 6.0. 
Both S and A are small at low pressures characteristic of the I phase. Here we adopt a value 
S ~ 0.4 as a heuristic threshold for the N phase Q, Q- From Fig. [3] it therefore appears 
that the N phase forms somewhere above P\\ ~ 1.0 (T = 4.0) and 1.2 (T = 6.0), respectively. 

At T = 4.0 there is a clear sign of a RN phase as S increases and A reaches a maximum 
value and then drops to small values for P\\ ~ 2.3; thus, a sequence of phases I-N-Sa-RN 
is found at T = 4.0 as P increases. At T = 6.0, A never rises above the residual value of 
about 0.1. Hence, the Sa phase does not form at this T. Nevertheless, we notice a small, 
step-like increase at P\\ ~ 2.2 in the plot of A. At this and all larger pressures considered 
the nematic order is high and increases even further reflected by a monotonic increase of S 
toward its limiting value 1.0. For P\\ > 2.2 the conflned fluid exhibits structural features of 

n 

the RN phase (see [131]). Three characteristic snapshots of the system are shown in Fig. HI 
The increasing orientational order of the system as P\\ increases is apparent. Incidentally, 
we note that the seemingly value 5* = 1 in the RN phase is a flnite-size effect, as was 
demonstrated by Eppenga and Frenkel 80 1. 



A question naturally arises: What is the relation between the RN phase found in this 
simulated conflned fluid and the bulk RNs? Our simulations show that the phase diagram 
of the conflned system is qualitative similar to phase diagram of the bulk fluid. The only 
change induced by the conflning surfaces is to lower the pressure of the phase transitions 
with respect to the bulk system. This shift in the phase diagram is a well-documented 
conflnement effect. Here, we consider spatial correlations between the centers of mass of 
molecules in the direction rj^ perpendicular to the orientation of a reference molecule, where 

= — {u-r)u. The perpendicular radial distribution function (y'(rx) is the probability to 
flnd a molecule at distance = \r±\ from a reference molecule. In Fig. [5]we compare g{r_i_) 
for a bulk and a conflned system in the RN phase. It is apparent that the maxima of g{r±) 
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FIG. 4: (Color online) "Snapshots" of characteristic configurations at T = 4.0; (a) Py = 1.6 (N), 
(b) P|| = 2.1 (Sa ), (c) P|| = 2.9 (RN) (see Fig. [3]). Molecules are aligned with n parallel to the x 
axis, whereas flat surfaces on top and bottom in dark represent the solid substrates. 

for the confined system are shifted to lower values of r_|_ with respect to the bulk case. This 
shift reflects the substrate-assisted efficient packing of molecules and therefore explains the 
shift of the Sa-RN transition to lower P\\ than in the bulk. However, the structure of the 
confined and bulk RN phases are very similar. 

Next, we address the important question of dynamics in the RN phase. To this end, 
we considered a specialized mean square displacement (MSD) to calculate displacements of 
molecules in the direction of their long axes. Specifically, we define r| = ■ and the 
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FIG. 5: (Color online) Perpendicular radial distribution function for a bulk (solid line) and confined 
(dashed line) simulation. Both states are in the RN phase at T = 4.0. 



associated MSD in the direction of the molecular long axes 

N 



\i=l 



r) - r|l (t) 



(15) 



From Eq. [15] a diffusion coefficient can be extracted using an Einstein relation, namely 

1 



Dh = lim — (Arf? (r)) 



(16) 



Figure [6] shows a dramatic increase of D\\ as one enters the RN phase whereas lower- 
pressure (I, N, or Sa) phases exhibit rather small self-diffusivity. In particular, the Sa phase 
is characterized by nearly vanishing self-diffusion constants which can be rationalized as 
above where we argued that the relatively compact layered structure makes it difficult for 
molecules to diffuse out of their original layer and penetrate into a neighboring one. The 
dramatic increase in mass transport in the direction of n in combination with nearly perfect 
nematic order prompted us to refer to liquid crystals in the RN phase as "supernematics" 



121. 



There are, in fact, experimental observations that are consistent with our simulations. 
For example, by extracting longitudinal relaxation rates from data of Miyajima et al. 32| . 
one finds that in the RN phase these relaxation rates are considerably lower than those char- 
acteristic of the N phase. At the same T, relaxation rates can be converted into correlation 
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FIG. 6: (Color online) Parallel self-diffusion coefficient as a function of transverse pressure Py. 
(a) Calculations at T = 4.0; (b) Calculations at T = 6.0. 

times that are significantly shorter in the RN as opposed to the N phase. However, this re- 
mains speculative, until experiments prove that translation diffusion is the main relaxation 
mechanism at this T. Stronger evidence is found in the conductivity experiments performed 



by Ratna [l6| (see Sec. IIIip . 



One may ask: Why is the diffusivity so large in the direction of the long molecular 
axis? As interesting analogy jl^ can be drawn with the case of diffusion in systems such as 
zeolites or single-wall carbon nanotubes. When the diameter of a diffusing particle matches 
the size of the confining nanochannel then the diffusivity increases markedly. This effect is 
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known as "levitation" Il5| . In our case, the high degree of orientational order in the 
RN phase causes a molecule to be trapped within the "cage" formed by its first-neighbors 
which, because of the prolate geometry, is effectively a narrow channel. Furthermore, the 
attractive intermolecular forces are overcome because the fluid is in a dense state (either low 
T or high pressure). Thus, within this channel a molecule fits the condition of the levitation 
effect. 



VI. CONCLUSIONS 

In this work, we have reviewed experimental and theoretical investigations of structural 
and dynamical properties of RN phases. It is surprising that after almost 40 years of research 
many details of RN phases are still escaping a comprehensive understanding. The main 
difficulty is the variety of substances exhibiting RN's. Pure substances and mixtures alike 
have been found with a RN in their phase diagrams. Polar and nonpolar molecules exhibit 
reentrancy alike. It is clear to the present authors that more experiments are necessary to 
settle some important, yet still open questions. For example: (i) How does the nematic 
order change upon crossing the Sa-RN transition? (ii) Is frustration really the mechanism 
underlying the reentrant behavior in LC's with a bilayer S phase? (iii) Experimentally, what 
is the mechanism driving reentrancy in LC's that do not form dimers? (iv) Do RN phases of 
single- component, nonpolar molecules really exhibit the huge self-diffusivity in the direction 
of n predicted by our simulations Q? 

Regarding the dynamical behavior, it is at present safe to assume that though many LC 
systems have a RN phase, their dynamics need not be alike in any respect. The reason 
is that different microscopic mechanisms may be responsible for reentrancy in different LC 
systems depending on the details of the intermolecular interaction potential. At the time of 



writing some experimental evidence exists 
a positive answer. Moreover, our findings 
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12 



suggesting that last question may indeed have 



13| seem consistent with an extrapolation of 



longitudinal relaxation rates of NMR experiments 32 1 . 
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